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Dynamics of the Density of Quantized Vortex-Lines in Superfluid Turbulence 
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The quantization of vortex lines in superfluids requires the introduction of their density C{r, t) in 
the description of quantum turbulence. The space homogeneous balance equation for C{t), proposed 
by Vinen on the basis of dimensional and physical considerations, allows a number of competing 
forms for the production term V. Attempts to choose the correct one on the basis of time-dependent 
homogeneous experiments ended inconclusively. To overcome this difficulty we announce here an 
approach that employs an inhomogeneous channel flow which is excellently suitable to distinguish 
the implications of the various possible forms of the desired equation. We demonstrate that the 
originally selected form which was extensively used in the literature is in strong contradiction with 
our data. We therefore present a new inhomogeneous equation for C{r,t) that is in agreement with 
our data and propose that it should be considered for further studies of superfluid turbulence. 


Background: Below the Bose-Einstein condensation 
temperature T\ Ri 2.18 K, liquid '^He becomes a quantum 
inviscid superfluidi^—. Aside from the lack of viscosity, 
the vorticity in ^He is constrained to vortex-line singu¬ 
larities of fixed circulation k = h/M, where h is Planck’s 
constant and M is the mass of the ^He atom. These 
vortex lines have a core radius oq ^ 10“® cm, compati¬ 
ble with the interatomic distance. In generic turbulent 
states, these vortex lines appear as a complex tangle with 
a typical intervortex distance . 

Recent progress in laboratory^^— and numerical stud¬ 
ies— of superfluid turbulence in superfluid ^He and 
“^He led to a growing consensus that the statistical prop¬ 
erties of superfluid turbulence at large scales R ^ £ are 
similar to those of classical turbulence. An acceptable 
theory of these large scale properties— is based on 
the Landau-Tisza two-fluid model— using ‘normal’ 
and ‘superfluid’ components of densities pn and ps with 
velocity fields Un{r,t) and Us{r,t). This model was ex¬ 
tended by Hall-Vinen^ and Bekarevich-Khalatnikov^ to 
include a mutual friction between the components, pro¬ 
portional to the the vortex-line density C. This means 
that a theory of large-scale motions which is affected by 
the mutual friction requires the inclusion of the dynamics 
of £. 

The situation changes drastically upon considering the 
statistical properties on smaller scales, where the quanti¬ 
zation of vortex lines becomes crucial. Several statistical 
characteristics of the vortex tangle become essential for 
a consistent description. Besides £, these characteris¬ 
tics involve the mean-square curvature (S')^, the vortex 
tangle anisotropy parameters /|| and R introduced by 
Schwara^. The most important of these is the vortex line 
density C{r,t). It is expected that coupling one or more 
of these quantities to the variables of the two fluid Hall- 
Vinen-Bekarevich-Khalatnikov equations—!^ is a mini¬ 
mal requirement for an acceptable theory of quantum tur¬ 
bulence. 

The problem: A phenomenological equation of mo¬ 
tion for C{t) was suggested by Vinen for homogeneous 


counterflows^Sii^ already in 1957: 

dC{f)/dt = 'P{f) — Vit) . (la) 

Here the production term P(t) describes the growth of C 
due to the extension of the vortex rings by mutual friction 
which is caused by the difference between the velocities of 
the normal and super components (“the counterflow ve¬ 
locity” 14s)- The decay term, 'D{f) is again caused by the 
mutual friction due to the moving normal fluid compo¬ 
nents and is assumed to be independent of 14s- Therefore 
both terms should be proportional to the dimensionless 
dissipative mutual friction parameter a. In principle it 
is not guaranteed that the equation for L can be closed 
via C and 14s- Such a closure for V and assumes 
that (beside a, k and 14s) the only relevant variable in 
the problem is the instantaneous value L{f), while S, /||, 
li, etc. are to some extent unimportant. Upon accepting 
this closure idea the dimensional reasoning dictates^ : 

r^Vci = a kC^F(x) , X = V^Jk^C , (lb) 
V=^V^i = anC^Gix) . (Ic) 

Here F{x) and G{x) are dimensionless functions of the 
dimensionless argument x. 

The most delicate issue in this approach is the determi¬ 
nation of the functions F{x) and G{x). Vinen^ assumed 
that the decay term D is independent of Vns leading to 
G{x) = Cdec- On the other hand Vinen and later authors 
(see, e.g— ) chose V to be proportional to the mutual 
friction force / oc a|14s|j leading to the proposition that 
F(x) oc -yT, and then 

Vci=^Vi=aGiC^/^\V^,\, (2a) 

where Gi is a dimensionless constant. Viner4^ realized 
that (l2a)l is not the only possibility. Another choice can 
follow the spirit of Landau’s theory of critical phenom¬ 
ena, considering C as an order parameter which deter¬ 
mines dC/dt via an analytical function F{x). Then the 
leading term in the expansion of F(x) oc x giving: 

Vci^V2 = aC2CV^jK . 


(2b) 
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FIG. 1; Numerical setup and obtained profiles of counterflow velocity and vortex line density. Panel a: Upper figure: Plane 
channel flow geometry Lx x Ly x L^, = 4:h,Ly = 2h,Lz = 2h,h = 0.05 cm. Lower figure: Projection of a vortex tangle 

on the [x, y)-plane, orthogonal to the walls. Different colors are used to distinguish vortex lines. Panels b and c: Prescribed 

parabolic (b) and non-parabolic (c) normal velocity profiles Un (j/^) (— • —), the resulting counterfiow profiles (-) and 

the resulting profiles of >C^(2/^) (—) in dimensionless unites defined by Eqs. ([8]). T = 1.6 K. 


Both options ([2]) are of course dimensionally correct and 
they predict the same stationary solution, £st oc 1/2, 
which is well supported by both experiments and nu¬ 
merical simulations (see, e.gj^ and references therein). 
In principle, one could hope to distinguish between the 
different forms of this important equation by compar¬ 
ing their prediction for the time evolution from some 
initial condition toward /Igt in the presence of counter- 
flow Vns. Unfortunately, the difference in prediction is 
too small. Neither Vinen^ himself nor later— experi¬ 
mental attempts succeeded to distinguish between these 
two forms—. We have made our own attempts to distin¬ 
guish between the two discussed models by numerical 
simulation of space homogeneous counterflow turbulence, 
finding inconclusive results as well. 

The proposed resolution of this old conundrum can 
be obtained by studying inhomogeneous flows like chan¬ 
nel flows in which the various relevant variables have non¬ 
trivial profiles. We will argue that in fact none of the 
equations m are correct. We propose yet a third form 
of V [corresponding to f{x) oc 

Pci =^V 3 = a Cp.odVcvi/K^ . (3a) 

Being dimensionally correct this closure fits the data that 
are presented below significantly better than either of the 
equations @. We are led to a revision of the homoge¬ 
neous equation of motion for the field C{r, f) in the form 

+ V • Jc\{r, t) = Vair, t) - Vcfir, f) . (3b) 

ot 

Here we have added a vortex-line density flux J'frfi). 
Based on our numerical simulations (see below) we sug¬ 
gest to model J'{r,t) as follows: 

J,fir,t) = -Cnn.{a/2n)VV^, . (3c) 

Notice that the suggested Eqs. m are based on our 
analysis of counterflow turbulence with laminar nor¬ 
mal fluid components. Nevertheless we believe that 


Eqs. m or their modifications may serve as a basis for 
future studies of inhomogeneous superfluid turbulence in 
a wide variety of conditions, including the evolution of a 
neutron-initiated micro big bang in superfluid ^He^, tur¬ 
bulent counterflows and pressure-driven superfluid chan¬ 
nel and pipe flows o/'^He, turbulent flow o/^He in rotat¬ 
ing cryostat^, etc.. 

Starting from first principles-. To reach the de¬ 
sired equation (l3b)) we denote the coordinates of the 
quantized vortex lines by s{^,f), which is parameterized 
by the arc-length Schwarz— derived the equation of 
motion for the length of the vortex-line segment <5^: 


1 dSi 
dt 


aVnsis,t) ■ {s' X s") . 


(4a) 


Here s' = ds/df,, s" = dfs/d^"^, a is the temperature 
dependent dissipative mutual friction parameter. The 
counterflow velocity Vns (which is a function of s and t 
which we suppress for notational simplicity) is the dif¬ 
ference between the normal fluid velocity V" and super¬ 
fluid velocity V®: 


Ks = V" - 1^% -V** = + Vbs . (4b) 


The super-fluid velocity includes the macroscopic po¬ 
tential part and the Biot-Savart velocity V^s- The 
later term is defined by the entire vortex tangle configu¬ 
ration C : 


Vb. 


K /" (s — Si) X dsi 

47r ic 


= V" 

r T T 


V; 


nl • 


(4c) 


The integral is logarithmically divergent when Si —>■ 
s. It is customary to regularize it by using the vor¬ 
tex core radius jio and the mean vortex line curvature 
radius R = 1/S. The main contribution to Vbs orig¬ 
inates from integrating over scales between oq and R, 
i.e. tto ^ jsi — s\ ^ R. This contribution is known as the 















3 





FIG. 2: For parabolic normal velocity profile: A comparison between production {V), decay (H) and flux {J) as obtained 
from numerical simulation and those obtained using their closure versions. Comparison of the numerical data ( —) to the 
competing production forms of the Vinen equation (Panel a): (— • —) -Vi, Eg. (I2al) . (—) - V 2 , Eg. (I2bl) and (— —) - P 3 , 
Eq. (I3bl) . Comparison of the numerical data for decay (Panel b) and flux term (panel c) with Eq. (IlcD for T> and Eq. (l3cD for J 
-(- -)• 


“Local Induction Approximation” (LIA)^ and is written 
as: 


= Ps' X s" , /3 = («/4^) In (i?/ao) . (4d) 

The term is non-local, being produced by the rest of 
the vortex configuration, C', with |si — s| > R -. 


K f (s — Si) X dsi 

4 tt Jc |s-SiP 


(4e) 


The next step is to integrate Eq. (I4al) over the vortex 
tangle in a fixed volume fl which resides in slices between 
y and y -I- 6y, going over all x and z. This provides us 
with the equation of motion for £{y, t) = This 

equation is written in the form similar to Eq. (I3bl) : 


dC{y,t) dj{y,t) 
dt dy 


'P{y,t) -V{y,t), 


(5) 


with the following identification for the flux J (toward 
the walls), production V and decay term V: 


Viy, t) = ^ f de (F" - - K^i) • X s"), (6a) 

“ JCn 

^ f = apcs'^ , (6b) 

“ JCn 

^ [ d^Vdrih,y = ^ [ dCKs.x4- (6c) 

“ -ICsi 

Here the production and decay terms come directly from 
integrating Eq. (I4al) ; they coincide with the correspond¬ 
ing equations in Ref.— with the only difference that in 
Eqs. (|n]) the integrals are taken in the slice il (between y 
and y -I- 5y)^ while in^ the integrals are taken over the 
entire volume. In the flux term the drift velocity Vdi.ift(^) 
of the vortex line segment s(^),, can be found from the 
vortex filament equations^ and written as follows: 


Kirift(0 = + as' X K, . (7) 


The mean value of V® is oriented along the x direction 
and it does not contribute to the y component of the flux, 
J. The j/-component of the second term in Eq. ([71) gives 
the final expression in Eq. (I6cll . 

Numerical simulations were set up in a 3- 
dimensional planar channel geometry (see Fig.|l}i) of half¬ 
width h with prescribed time-independent profile of the 
streamwise projection of the normal velocity V^{y). To 
find the vortex tangle configurations we used the vortex 
filament method, taking into account the potential flow 
to maintain the counterflow condition. Details of the 
simulation method can be found in Refs.— Here we 
used the reconnection method^ and the line resolution 
= 1.6 X 10“^ cm. Periodic conditions were used in the 
streamwise and spanwise directions. Taking into account 
the fact that the boundary conditions for the superfluid 
component are still under discussion we adopt their sim¬ 
plest version: in the wall-normal y direction Vy{±h) = 0 
and s'{zLh) = (0,±1,0) at the solid walls. Periodically 
wrapped replicas of the tangle configuration were used in 
the X- and z- directions, with reflected configurations in 
the y direction. 

Having selected a stationary profile of V^{y) = Vn{y) 
we started with a set of arbitrary oriented circular vor¬ 
tex rings and solved the equation for the vortex line evo¬ 
lution. For the obtained dense vortex tangle we found 
the time-averaged profile {V^{y))^ = Vs{y). In all cases 
we ran the simulations until we obtained steady mean 
profiles. We take temperature T = 1.6 K with mutual 
friction coefficients a = 0.098, a' = 0.016, Ref—. 

We begin with a parabolic profile for Vn{y). The pro¬ 
files of Vns{y) and C{y) in the dimensionless form; 

yt = y/L , Rt = E/v^, £1 = ^c/ ) , (8) 

are shown in Fig.|T[i. Taking integrals in Eqs. ([5]) over 
the numerically found vortex tangle configuration we can 
compute the production, decay and flux terms, denoted 
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FIG. 3: For non-parabolic normal velocity profile: A comparison between production {V), decay {V) and flux {J) as obtained 
from numerical simulation and those obtained using their closure versions. Comparison of the numerical data ( —) to the 
various allowed production forms of the Vinen equation (Panel a): (— • —) -Pi, Eq. (|2a|l . (—) - P 2 , Eq. (I2b|l and (— —) - P 3 , 
Eq. (I3bl) . Comparison of the numerical data for decay (Panel b) and flux term (panel c) with Eq. (Ilcl) for D and Eq. (I3cl) for J 
-(- -)• 


as 'Pnum{y,t), 'Dnuni{y,t) and Jnum{y,t), respectively . 
Then we compared them with the various their closure 
versions, Vi, V 2 , Vch and Jci- 

From the theoretical point of view the questions are: 
can we approximate the integrals in Eqs. (jH]) only in terms 
of the counterflow velocity 14s and C itself, or would the 
integrals produce other dynamical variables that should 
require further coupled equations to close the system? Is 
closure possible in general, or only in some conditions? 

Assuming that closure is allowed, dimensional consid¬ 
erations presented us with different versions for the pro¬ 
duction term Vn{y) given by Eqs. (|2al) . (iTbll and (I3a1) for 
n = 1,2 and 3. With the widely accepted approxima¬ 
tion— that the tangle curvature radius i? = I/S' is pro¬ 
portional to the intervortex distance £ = 1 /'/L one gets 
from Eqs. (IFbl) and (lidl) the closure (ITB for V with 

Cdec = C 2 ln(i?/ao)/47r , where C 2 = iS. (9) 

For the flux term ((5^ we suggest (in the channel geome¬ 
try): 

Jci{y,t) = -{a/2K)CfiuKdV^Jdy . (10) 

In Fig. [2] we compare the numerical integrals® (shown 
as thick solid black lines) with corresponding closures. 
The dot-dashed blue line in Fig. [5^. shows the Vinen pre¬ 
diction, Pi(?/), (I2a1) . while the thin solid green line shows 
the alternative form P 2 (y), Eq- (l2bl) . By dashed red line 
we show the prediction P 3 (y),which is evidently superior 
to the other two. From this data we can conclude that 
Eq. (iTbl) is the one that should be used in the present 
inhomogeneous case. 

Figure[2]3 shows that the numerical integral (l6al) and 
the commonly used form (O for the decay term, V, prac¬ 
tically coincide, meaning that C 2 , defined by Eq. ®, is 
indeed ^-independent. Figure [SJ: also demonstrates very 
good agreement between Jnumiy) and Jci{y) given by 
Eqs. (15^ and (ITUl) . 


K ^/WS/y^ C-p.od Cdec Gflux Cdec (V/s) Cflux/ {Vl) 
cm/s (cm/s)^ (s/cm)^ 


1 

1.0 

0.788 

0.018 

6.9 

0.088 

4.31 

0.145 

2 

1.2 

0.785 

0.018 

4.9 

0.124 

4.35 

0.147 

3 

1.5 

0.780 

0.019 

3.2 

0.210 

4.32 

0.146 

4 

1 

1.247 

0.022 

2.3 

0.048 

3.62 

0.039 


TABLE I: Closure constants for different values of the nor¬ 
mal velocity maximum 14 for parabolic (lines 1-3) and non¬ 
parabolic (line 4) profiles. 


Realizing that the good match between the numerical 
data and Eq. (I3bl) may be accidental due to particulary 
chosen numerical parameters we repeated the simulations 
with other magnitudes of the counterflow velocities. We 
found again a good agreement between the numerical 
profiles Pnum(y), Vnumiy) and Jnum(y) with the corre¬ 
sponding closures, Vsiy), 'Ddiy) and Jci{y)- Consistency 
requires that the numerical constants Cprod, Cdec and 
Cflux in these closures should be 14s independent. Table 
Ushows that this is the case only for Cprod, while Cdec and 
Cflux approximately depend on 14s as: Cdec oc 1/ (Ens) 
and Cflux oc (Eng), where (l4^s) is the mean-square of 
14s across the channel (which, in its turn oc E,/). From 
these facts we can conclude that our closure (153 for the 
production term is confirmed, while the traditional clo¬ 
sure © for the decay term and simple closure (HOD seems 
to be questionable, although they reproduce well the pro¬ 
files Vnumiy) and J^numiy) in the parabolic case. 

To clarify further the quality of the discussed closures 
we imposed non-parabolic normal velocity profile with 
two maxima and zero on the centerline, shown in Fig.[T}:. 
Although this profile looks strange, it may be realized in 
a counterflow experiment with non-homogeneous heat¬ 
ing in short enough channel^. For this profile we again 
computed the counterflow and vortex-line density pro- 
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files shown in Fig.[T]:. Next, we repeated all the steps 
described before, and found again that our proposed 
form Eq. (|3b| fits the data much better than the other 
two forms as seen in Fig. [5^. Figures[5}D,c show that the 
closures for the decay and the flux terms, 2?ci and Jd, 
are running into trouble reflecting the numerical pro¬ 
files only very roughly. In particular this means that 
c| = which is involved in the closure (IT^ via 

Eq. ® , varies across the channel by a factor of three, as 
follows from Fig.|3jD. Recall, that in the parabolic case c| 
is y-independent, although it depends on (Fn^s) approxi¬ 
mately as 1/ (l4^s)- Consequently, to improve the closure 
for the decay term one needs to involve an additional bal¬ 
ance equation for the tangle curvature S. Similarly, our 
analysis shows that to improve the closure for the flux 
term, one needs to involve information about the tangle 


anisotropy. 

Conclusion: The suggested closure (I5al) for the pro¬ 
duction Pprod(2/) using the counterflow velocity and the 
vortex line density profiles can be considered as highly 
promising. However, the closures for the decay and flux 
are sensitive, and they may require accounting for addi¬ 
tional tangle characteristics. The first candidate is the 
tangle curvature; the tangle anisotropy also can be im¬ 
portant. Much more work in this direction is required to 
develop a consistent theory of the wall-bounded super¬ 
fluid turbulence. 
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knowledges the kind hospitality at the Weizmann Insti¬ 
tute of Science during the main part of the project. 


^ R. J. Donnelly, Quantized Vortices in Hellium II (Cam¬ 
bridge Univ. Press, Cambridge, UK, 1991). 

^ Quantized Vortex Dynamics and Superfluid Turbulence, 
ed. by C. F. Barenghi et al., Lecture Notes in Physics 571 
(Springer-Verlag, Berlin, 2001). 

® A. F. Annett, Superconductivity, Superfluids and Conden¬ 
sates (Oxford University Press, Oxford, 2004). 

^ L. Skrbek, K. R. Sreenivasan, Developed quantum turbu¬ 
lence and its decay, Phys. Fluids 24 011301(2012). 

® S. K. Nemirovskii, Quantum turbulence: Theoretical and 
numerical problems, Phys. Rep. 524, 85 (2013). 

® C. F. Barenghi, V. S. L’vov and P.-E. Roche, Experimen¬ 
tal, numerical, and analytical velocity spectra in turbulent 
quantum fluid, Proc. Natl. Acad. Sci. USA 111, (Supl.l), 
4683 (2014). 

^ C. F. Barenghi, L. Skrbek and K. R. Sreenivasan, Introduc¬ 
tion to quantum turbulence, Proc. Natl. Acad. Sci. USA 
111, (Supl.l), 4647 (2014). 

® V.B. Eltsov, R. de Graaf, R. Hanninen, M. Krusius, R.E. 
Solntsev, V.S. L’vov, A.I. Golov, P.M. Walmsley, Turbu¬ 
lent dynamics in rotating helium superfluids, Progress in 
Low Temperature Physics, XVI pp. 46-146 (2009). 

® E. Fonda, D. P. Meichle, N. T. Ouellette, S. Hormoz, and 
D. P. Lathrop, Direct observation of Kelvin waves excited 
by quantized vortex reconnection, Proc. Natl. Acad. Sci. 
USA 111 , (Supl.l), 4707 (2014). 

W. Guo, M. La Mantiac, D. P. Lathrop and S.W. 
Van Sciver, Visualization of two-fluid flows of superfluid 
helium-4, Proc. Natl. Acad. Sci. USA 111, (Supl.l), 4653 
(2014). 

L. Skrbek, A.V. Gordeev, F. Soukup, Decay of counter- 
flow He II turbulence in a finite channel: Possibility of 
missing links between classical and quantum turbulence, 
Phys. Rev. E 67, 047302 (2003). 

S. W. Van Sciver, S. Fuzier, and T. Xu, J., Particle Im¬ 
age Velocimetry Studies of Counterflow Heat Transport in 
Superfluid Helium H, J. Low Temp. Phys. 148, 225 (2007). 

T. Zhang and S. W. Van Sciver, Large-scale turbulent flow 
around a cylinder in counterflow superfluid 4He (He (II)), 
Nat. Phys. 1, 36 (2005). 

G. P. Bewley, D. P. Lathrop, and K. R. Sreenivasan, Super- 
fluid Helium: Visualization of quantized vortices. Nature 


(London) 441, 588 (2006). 

M. S. Paoletti, R. B. Fiorito, K. R. Sreenivasan, and D. P. 
Lathrop,Visualization of Superfluid Helium Flow, J. Phys. 
Soc. Jap. 77, 111007 (2008). 

D. N. McKinsey, W. H. Lippincott, J. A. Nikkei, and 
W. G. Rellergert, Trace Detection of Metastable Helium 
Molecules in Superfluid Helium by Laser-Induced Fluores¬ 
cence, Phys. Rev. Lett. 95, 111101 (2005). 

W. G. Rellergert, S. B. Gahn, A. Garvan, J. G. Hanson, W. 
H. Lippincott, J. A. Nikkei, and D. N. McKinsey, Detection 
and Imaging of He2 Molecules in Superfluid Helium, Phys. 
Rev. Lett. 100, 025301 (2008). 

W. Guo, J. D. Wright, S. B. Gahn, J. A. Nikkei, and D. 

N. McKinsey, Metastable Helium Molecules as Tracers in 
Superfluid He4, Phys. Rev. Lett. 102, 235301 (2009). 

W. Guo, S. B. Gahn, J. A. Nikkei, W. F. Vinen, and D. N. 
McKinsey, Visualization Study of Gounterflow in Super- 
fluid He4 using Metastable Helium Molecules Phys. Rev. 
Lett. 105, 045301 (2010). 

A. Marakov, J. Gao, W. Guo, S. W. Van Sciver, G. G. 
Ihas, D. N. McKinsey, and W. F. Vinen, Visualization of 
the normal-fluid turbulence in counterfiowing superfluid 
He4. Phys. Rev. B 91, 094503 (2015). 

G.F. Barenghi, A.V. Gordeev, L. Skrbek, Depolarization 
of decaying counterflow turbulence in He H, Phys. Rev. E 
74, 026309 (2006). 

M. Sciacca, Y.A. Sergeev, C.F. Barenghi, L. Skrbek, Sat¬ 
uration of decaying counterflow turbulence in helium H, 
Phys. Rev. B 82, 134531 (2010). 

N. G. Berloff, M. Brachet and N. P. Proukakis, Modeling 
quantum fluid dynamics at nonzero temperatures, Proc. 
Natl. Acad. Sci. USA 111, (Supl.l), 4675 (2014). 

R. Hanninen and A. W. Baggaley, Vortex filament method 
as a tool for computational visualization of quantum tur¬ 
bulence, Proc. Natl Acad. Sci. USA 111, (Supl.l), 4667 
(2014). 

A. W. Baggaley and Laizet, Vortex line density in coun¬ 
terfiowing He II with laminar and turbulent normal fluid 
velocity profiles, Phys. Fluids 25, 115101 (2013). 

L. Kondaurova, V. L’vov, A. Pomyalov and 1. Procaccia, 
Structure of a quantum vortex tangle in ^He counterfiow 
turbulence, Phys. Rev. B 89, 014502 (2014). 



6 


G. V. Kolmakov, P. V. E. McClintock and S. V. Nazarenko, 
Wave turbulence in quantum fluids, Proc. Natl. Acad. Sci. 
USA 111 , (Supl.l), 4727 (2014). 

H. Adachi, S. Fujiyama, M. Tsubota, Steady-state counter¬ 
flow quantum turbulence: Simulation of vortex filaments 
using the full Biot-Savart law, Phys. Rev. B 81, 104511 
( 2010 ). 

V. S. L’vov, S. V. Nazarenko and G. E. Volovik, Energy 
spectra of developed superfluid turbulence, J. Low Temp. 
Phys. 80, 535 (2004). 

V.S. L’vov, S.V. Nazarenko and L. Skrbek, Energy Spectra 
of Developed Turbulence in Helium Superfluids, J. Low 
Temp. Phys. 145, 125 (2006). 

V. S. L’vov, S. V. Nazarenko and O. Rudenko, Gradual 
eddy-wave crossover in Superfluid turbulence, J. of Low 
Temp. Phys. 153, 140-161 (2008). 

L. Boue, V.S. L’vov, A. Pomyalov, 1. Procaccia, Energy 
spectra of superfluid turbulence in He-3B, Phys. Rev. B 
85, 104502 (2012). 

L. Boue, V.S. L’vov, A. Pomyalov, 1. Procaccia, Enhance¬ 
ment of intermittency in superfluid turbulence, Phys. Rev. 
Lett., 110, 014502 (2013). 

L. Boue, V S. L’vov, Y. Nagar, S. V. Nazarenko, A. 
Pomyalov, and 1. Procaccia, Energy and Vorticity Spec¬ 
tra in Turbulent Superfluid 4He from T = 0 to Ta, Phys. 
Rev. B 91, 144501 (2015). 

L.D. Landau, Theory of superfluidity of helium-II, J. Phys. 
USSR 5, 71 (1941). 

L. Tisza, J. Phys. Radium 1, 164 (1940) , ibid 1, 350 
(1940). 

H.E. Hall, W.F. Vinen, The rotation of liquid helium H. H. 


The theory of mutual friction in uniformly rotating helium 

H, Proc. R. Soc. Lond. A Math Phys. Sci. 238, 215 (1956). 

I. L. Bekarevich, I.M. Khalatnikov, Phenomenological 
Derivation of the Equations of Vortex Motion in He H, 
Zh. Eksp. Teor. Fiz. 40, 920 (1961) (Sov. Phys. JETP 13, 
643 (1961)). 

K.W. Schwarz, Three-dimensional vortex dynamics in su¬ 
perfluid ^He: Homogeneous superfluid turbulence, Phys. 
Rev. B 38, 2398 (1988). 

W. F. Vinen, Mutual friction in a heat current in liquid 
helium H. 1. Experiments on steady heat currents, Proc. 

R. Soc. Lond. A Math. Phys. Sci. 240, 114 (1957). 

N. B. Kopnin. Vortex Instability and the Onset of Super¬ 
fluid Turbulence Phys. Rev. Lett. 92, 135301 (2004). 

W. F. Vinen, Mutual friction in a heat current in liquid he¬ 
lium 11. 11. Experiments on transient, Proc. R. Soc. Lond. 
A Math. Phys. Sci. 240, 128 (1957). 

S. K. Nemirovskii and W. Fiszdon, Ghaotic quantized vor¬ 
tices and hydrodynamic processes in superfluid helium. 
Rev. Mod. Phys. 67, 37 (1995). 

Y.M. Bunkov, A. 1. Golov, V. S. L’vov, A. Pomyalov and 1. 
Procaccia, Evolution of Neutron-Initiated Micro-Big-Bang 
in superfluid He 3B, Phys. Rev B 90, 024508 (2014). 

D. G. Samuels , Velocity matching and Poiseuille pipe flow 
of superfluid helium, Phys. Rev. B 46, 11714 (1992). 

R. J. Donnelly , C. F. Barenghi , The Observed Properties 
of Liquid Helium at the Saturated Vapor Pressure, J. Phys. 
Ghem. Ref. Data 27, 1217(1998). 

Skrbek, L., private communication. 

For more detailed discussion of this problem see, e.g.— 


